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Abstract 

The states of an open quantum system are coupled via the environment of scattering wave- 
functions. The complex coupling coefficients uj between system and environment arise from the 
principal value integral and the residuum. At high level density where the resonance states overlap, 
the dynamics of the system is determined by exceptional points. At these points, the eigenvalues 
of two states are equal and the corresponding eigenfunctions are linearly dependent. It is shown in 
the present paper that Im(cj) and Re(w) influence the system properties differently in the surround- 
ing of exceptional points. Controlling the system by a parameter, the eigenvalues avoid crossing 
in energy near an exceptional point under the influence of Re(w) in a similar manner as it is well 
known from discrete states. Im(u;) however leads to width bifurcation and finally (when the system 
is coupled to one channel, i.e. to a common continuum of scattering wavefunctions), to a splitting 
of the system into two parts with different characteristic time scales. Physically, the system is 
stabilized by this splitting since the lifetimes of most (N — 1) states are longer than before while 
that of only one state is shorter. In the cross section the short-lived state appears as a background 
term in high-resolution experiments. The wavefunctions of the long-lived states are mixed in those 
of the original ones in a comparably large parameter range. Numerical results for the eigenvalues 
and eigenfunctions are shown for A = 2, 4 and 10 states coupled mostly to 1 channel. 
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I. INTRODUCTION 



The description of quantum mechanical systems by means of the Schrodinger equation 
has been developed more than 80 years ago. At that time only a few resonance states in 
nuclei and atoms were known which are well separated from one another. The energies 
of these states are well described by means of the Schrodinger equation with a Hermitian 
Hamiltonian. In order to describe the finite lifetimes of these states, the R-matrix theory 
has been developed which is however too complicated for practical calculations. The finite 
lifetimes of the individual states are calculated usually perturbatively, see e.g. jl]. Later, 
systems at high level density were in the center of interest. Under these conditions, a 
statistical description of the system turned out to be very efficient |2|. The resonances are 
well isolated from one another, and the lifetimes of the states are very long and do not play 
any role. 

In the course of time, experimental studies have been performed for different quantum 
systems with a much improved accuracy. Also the theoretical calculations are carried out 
today not only for states being well separated in energy from one another but also for 
resonance states in the regime of overlapping. Here the single individual states can no 
longer be identified what results in problems of their interpretation. Contradictions between 
experimental results and conventional Hermitian quantum physics appeared in different 
small quantum systems. For example, Heiblum et al. [3( found experimentally a crossover 
from the mesoscopic to a universal phase for electron transmission in quantum dots about 
than 10 years ago. These results could not be explained in conventional quantum physics 
in spite of much effort |4|. Recently Koehler et al. [5] observed that neutron resonance 
data exclude random matrix theory. Deviations between experimental data and random 
matrix theory in nuclear physics studies were observed also earlier, e.g. [6(. These and 
other experimental results show that the Schrodinger equation originally introduced with 
a Hermitian Hamilton operator for the description of well isolated resonances, has to be 
extended. Above all, the lifetimes of the resonance states have to be calculated also in the 
regime of overlapping resonances and the justification of statistical assumptions has to be 
proven for small quantum systems. 

The lifetimes of the resonance states can be calculated when the system is explicitly 
considered to be open and the calculations are performed quantum mechanically for both the 
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system and the environment of scattering wavefunctions into which the system is embedded. 

Q 

Using the Feshbach projection operator technique [7], first the energy-independent many- 
body problem of the system (with the Hermitian Hamiltonian Hb) is solved in the standard 
manner. These solutions provide the energies Ef and wavefunctions $f of the discrete 
states with inclusion of the so-called internal interaction. In a second step, the energy- 
dependent scattering wavefunctions £f of the environment are calculated and, further, the 
(energy-dependent) coupling matrix elements 

_ „ fnZ I^B\ir\cE 



iL = ^W\v\&). (1) 

between the discrete states of the system and the environment are evaluated (see BJ, section 
2.1). The corresponding Schrodinger equation contains an energy dependent (nonlinear) 
source term describing the coupling between system and environment. The Hamiltonian 
"Ho °f this equation is non-Hermitian and provides the lifetimes of the states, but without 
any additional interaction of the states via the environment. The Schrodinger equation 
with Hq and source term can be rewritten into a Schrodinger equation without source term 
but with a non-Hermitian Hamiltonian % that contains the interaction of the states via the 
environment (the so-called external interaction) in the nondiagonal matrix elements [9] . The 
Hamiltonian % reads 

U = H B + VbcG^Vcb (2) 

where Vbc an d Vcb stand for the coupling between system and environment and Gq is the 
Green function in the subspace of scattering states. The external interaction of the states 
via the continuum is complex, generally. The principal value integral is 

Re (4>f \U\®f ) - E?8 i3 = i- jy£ dE> (3 ) 
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and the residuum reads 

Im (#f \H\*f ) = ~ E7°7? c - (4) 

Z c=l 

The interaction of the states of the system via the environment is involved in the eigenval- 
ues Si and eigenf unctions of the Hamiltonian TL. It is therefore relatively easy, in this 
formalism, to study the influence of the environment onto the states of the system. That 
means, the non-Hermitian quantum physics is - in contrast to the widely spread meaning 



- not a further approximation introduced in the theory. It is an expression of the fact that 
the system considered is really open and its properties are influenced by the coupling to the 
environment. This influence is unimportant at low level density where it can be described by 
perturbation theory. It becomes, however, decisive in the regime of overlapping resonances. 

Meanwhile, many calculations for open quantum systems are performed with high accu- 
racy. Without using any perturbation theory or statistical assumptions, the Schrodinger 
equation with the non-Hermitian Hamilton operator % is solved and the eigenvalues 
Si = Ei — |Tj and eigenf unctions $j are obtained. That means, not only the energies 
Ei of the states of the system are evaluated but also their lifetimes T{ oc 1/IV The control of 
the eigenvalues by a parameter allows to draw conclusions on the dynamics of open quantum 
systems. 

The limits of the applicability of the standard quantum theory with Hermitian Hamilton 
operator can be seen best at high level density, where the individual resonance states overlap 
and interact with one another via the continuum of scattering wavefunctions. Here, the dy- 
namics of the system is determined by singular points, the so-called exceptional points. They 
cause level repulsion in energy as well as a bifurcation of the widths (inverse proportional 
to the lifetimes) of the states. Results are obtained theoretically as well as experimentally, 
which are counterintuitive at first glance : due to width bifurcation, coherent short-lived 
states are formed together with states that are almost decoupled from the continuum of 



scattering states 
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llj . An example is the above-mentioned crossover from the meso- 
scopic to a universal phase in the transmission in quantum dots 3j and its qualitative 
explanation on the basis of a Schrodinger equation with non-Hermitian Hamilton operator 



121 ] . This phenomenon is observed in different systems and is called usually dynamical phase 



transition [13j. The short-lived states are analog to the coherent superradiant states con- 
sidered by Dicke [14J] in 1954 as shown by means of a toy model 15]. Moreover, it could be 



shown in an experiment on non-locally coupled pairs of quantum point contacts that discrete 
states undergo a robust interaction that is achieved by coupling them to each other through 
the continuum Jl6]. Most of these results are of high value for fundamental questions of 
quantum mechanics as well as for applications. 

It is the aim of the present paper to show some generic (mostly numerical) results for 
open quantum systems in order to receive a deeper understanding for the dynamical phase 
transitions occurring in the regime of overlapping resonances. A toy model is used in order 
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to receive conclusions on the role played by exceptional points and avoided level crossings at 
high level density in open quantum systems. Of special interest are the effects arising from 
the imaginary part (jlj) of the coupling term via the environment. 

The calculations are performed with respectively two, four and ten resonance states cou- 
pled mostly to one open decay channel (corresponding to the common continuum of scat- 
tering wavefunctions) under the condition that they avoid crossing or even cross in a certain 
parameter range. The coupling of N resonance states to K < N channels corresponds to 
the general situation of open quantum systems. The formalism used in the calculations is 
the same as that discussed in our earlier paper [17]. The results show very clearly that 
the imaginary part of the coupling term between system and environment causes width 
bifurcation and, finally, the formation of different time scales in the regime of overlapping 
resonances. Due to width bifurcation, the system splits into two parts that exist at different 
times. This process is irreversible. The wavefunctions of the short-lived aligned states are 
mixed coherently in relation to the open decay channel while those of the long-lived trapped 
states are mixed incoherently such that the states are almost decoupled from the open decay 
channel. 

The results confirm the characteristic features of the dynamical phase transitions appear- 
ing in open quantum systems at high level density. They show that the number of states 
of the system is reduced during the dynamical phase transition due to ejecting the aligned 
short-lived state and, furthermore, that the wavefunctions of the states before and beyond 
the dynamical phase transition are completely different from one another. While the states 
of the original system (without interaction via the continuum) have individual spectroscopic 
features, those of the system consisting of the long-lived states beyond the dynamical phase 
transition show chaotic features. These results prove once more that the experimental re- 
sults [3j on the crossover from the mesoscopic to a universal phase for electron transmission 
in quantum dots can be explained by means of the non-Hermiticity of the Hamiltonian. 
Additionally it should be mentioned here that the aligned state corresponds to the Dicke 
superradiant state and the trapped states to the subradiant states which are considered in 



optics, e.g. 



18|- 



In section [Til we sketch the formalism used in the present calculations. For N = 2 
states with equal decay widths, the analytical and numerical solutions of the problem show 
clearly the effect of width bifurcation. Numerical results for eigenvalues and eigenfunctions 
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obtained for N = 2, 4 and 10 states are provided in sections II III to IVIl The results are 
discussed in section IVHI while conclusions on the relation between exceptional points and 
width bifurcation are drawn in the last section IVIIII 



II. FORMALISM 



We consider an N x N matrix 



ei + wn 



Wis 



\ 



H 



w 2 i 



62 + C022 



^>2N 



(5) 



UN2 




the diagonal elements of which are the N complex eigenvalues 6j + Ua = — i/2 ji of a 
non-Hermitian operator. The Ua are the so-called selfenergies of the states arising from their 
coupling to the environment of scattering wavefunctions into which the system is embedded. 
In atomic physics, these values are known as Lamb shift. Our calculations are performed with 
coupling matrix elements w„ the values of which do not depend on the parameter considered. 
In such a case, the Ua can considered to be included into the diagonal matrix elements, which 
read Si = ei + Uu. The e, and % denote the energies and widths, respectively, of the N states 
(including their selfenergies) without account of the interaction of the different states via 
the environment. 

The internal interaction of the two states % and k ^ i (appearing in the closed system) as 
well as their external interaction (via the environment) are contained in the io^. The internal 
interaction can be caused only by some part of Ke(uik) while the external interaction contains 
complex Uik, see equations (j3J) and OH). Most interesting part of the external interaction is 
therefore Im(a;jfc)- It becomes important at high level density where the corresponding 
resonance states overlap. 

When the number N of states is equal to the number K of common channels (i.e. equal 
to the number K of different common continua of scattering wavefunctions) all the coupling 
matrix elements cu^ are different from zero and the matrix (J5J) is full. In the case with 
only one open decay channel K — 1, all Uik different from Ui k=K and i^i=K k are zero. An 
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example of iV = 4 states coupled to only the fourth channel is the following matrix 

'ex u 14 ^ 



H 



V 



(6) 



/ 



e 2 w 2 4 
e 3 u u 

Ul 41 U42 U43 64 

with £j = €i for i 7^ 4 and £4 = €4 + W44. 

The eigenvalues of H will be denoted by Si = E$ — z/2 I\ where i% and Tj stand for the 
energy and width, respectively, of the eigenstate i. The eigenfunctions of the non-Hermitian 
~H are biorthogonal (see sections 2.2 and 2.3 of |8[), 



(*.*l*;> 



5ij . 



It follows 



and 



= Re ; ^ = > 1 

1511 = 1(^1^1 > 0. 



(7) 



(8) 



(9) 



At an exceptional point, Ai — > 00 and \B\ 



00. The Si and $j contain global features 



that are caused by many-body forces induced by the coupling uj^ of the states % and k ^ i 
via the environment, see equations (J3J), (jlj) and the corresponding discussion in jlll. 



In the case N = 2, the two eigenvalues of "H are 



Si j — Eij 



£1 +£2 



±Z; Z=-^/(ei 



s 2 ) 2 + 4cu 2 . 



(10) 



According to this expression, two interacting discrete states (with 7^ = 0) avoid always 
crossing since u and si — £2 are real in this case. Resonance states with nonvanishing widths 
Ti repel each other in energy according to the value of Re(Z) while the widths bifurcate 
according to the value of lm(Z). The two states cross when Z = 0. This crossing point is 
an exceptional point according to the definition of Kato 20] . 

At the exceptional point Z = 0, the eigenfunctions of fl5]) of the two crossing states are 
linearly dependent from one another, 



$7 -> ± i $ 



2 1 



(11) 



according to analytical as well as numerical and experimental studies, see Appendix of 
11] and section 2.5 of |8J. That means, the wavefunction $1 of the state 1 jumps, at the 
exceptional point, via the wavefunction $! ± i$ 2 of a chiral state to ±i$ 2 [21|. From 
(TUD follows : 

(i) When two levels are distant from one another, their eigenf unctions are (al- 
most) orthogonal, (<&* k \®k) ~ (®k\®k) = A k « 1. 

(ii) When two levels cross at the exceptional point, their eigenf unctions are 
linearly dependent according to ( TTTT) and ($fc|$fc) = A k — > oo. 

These two relations show that the phases of the two eigenfunctions relative to one another 
change when the crossing point is approached. This can be expressed quantitatively by 
defining the phase rigidity r k of the eigenfunctions <E>fc, 

(gg|gj) _ x 

rk ~ " Afe • (12) 

It holds 1 > r k > 0. The non-rigidity r k of the phases of the eigenfunctions of % follows 
also from the fact that ($£|<E>fc) is a complex number (in difference to the norm ($> k \<&j.) 
which is a real number) such that the normalization condition (J7]) can be fulfilled only by 
the additional postulation Im($£|$fc) = (what corresponds to a rotation). 

When rf. < 1, an analytical expression for the eigenfunctions as a function of a certain 
control parameter can, generally, not be obtained. An exception is the special case 71 = 72 



for which Z = |y (e\ — e 2 ) 2 + 4a; 2 . In this case, the condition Z = can not be fulfilled if 
lo = x is real due to 

( ei -e 2 ) 2 + 4x 2 > 0. (13) 

The exceptional point can be found only by analytical continuation into the continuum jsl, Q 
and the two states avoid crossing. This is analogue to the avoided level crossings of discrete 
states. 

The condition Z = can be fulfilled however if oj = i x is imaginary, 

(ei -e 2 ) 2 -Ax 2 = -»■ e 1 -e 2 = ±2x, (14) 

and two exceptional points appear. It holds further 

(ei -e 2 ) 2 > 4x 2 ->■ Z G « (15) 
(ei -e 2 ) 2 < 4a; 2 ^ Z G 3 (16) 
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independent of the parameter dependence e«(a). In the first case, the eigenvalues Si = 

Ei — ijlYi differ from the original values £j = e; — i/2 7, by a contribution to the energies 

and in the second case by a contribution to the widths. The width bifurcation starts at one 

of the exceptional points and becomes maximum in the middle between the two exceptional 

points. This happens at the crossing point e\ = e 2 where AT/2 = |Ti/2 — T 2 /2| = 4x. 

Some years ago, the case N = 2 with = e«(a), fixed real u = ujyi = CO21, and different 

fixed values of 7i, including 7, = 0, is studied as a function of the parameter a in the 

n 

neighborhood of avoided and true crossings of the two levels j9j. The results for the N = 2 
case 9s] show further that the wavefunctions of the two states $1 and $ 2 are mixed in a 
finite range of the parameter a around the critical value a cr at which the two states avoid 
crossing. This holds true not only for resonance states but also for discrete states. 

The eigenf unctions <£>j of H can be represented in the set of basic wavefunctions <pi of 
the unperturbed matrix (corresponding to the case with vanishing coupling matrix elements 

*i = EMj- ( 17 ) 

Also the bij are normalized according to the biorthogonality relations of the wavefunctions 

In our calculations, the mixing coefficients of the wavefunctions of the two states due 
to their avoided crossing are not calculated. We simulate the fact that the two wavefunctions 



are mixed in a finite parameter range around the critical value of their avoid crossing {9] by 
assuming a Gaussian distribution 

w ¥j = u e - (e '-^ )2 (18) 



for the coupling coefficients. The results reproduce very well those discussed in (9] for 2 
levels and real coupling u. 

Of special interest is the situation at high level density where the ranges of avoided 
crossings, defined by (Fl8|) . of different levels overlap. Some generic results obtained with 2, 
4 and 10 resonance states will be presented in the following sections. 
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III. EIGENVALUES AND EIGENFUNCTIONS FOR N = 2 CROSSING LEVELS 



We start our calculations with the two-level case which is studied mostly in literature. 
We choose the matrix §5§ with = ej(a), different fixed values of 7^ and fixed u. The 
functional dependence of the energies over the parameter a is similar as in 19, Ol. However, 

n n 

the uj are real in |9J while they are mostly complex in the present paper as well as in |17| . 
According to (fTOl) to f[T6|) . the exceptional points appear at different values of the parameter 
a when the ratio Re(u) to Im(w) is varied, see also Figs. 1 to 3 in 

In Fig. [TJwe show the numerical results for the avoided level crossing phenomenon as a 
function of the parameter a beyond, at and below the exceptional point. In all cases, the 
value of to is fixed to 0.05 (1 + i) while the 7i/2 are different in the different subfigures 
(with 72/71 = 1.1 in all cases). The exceptional point appears at 7i r /2 = 1, see Figs. [TJ c, 
d. When 7i/2 > 7^ r /2 (Figs. [TJ a, b), the eigenvalue trajectories cross in energy while the 
widths Ti and T 2 are always different from one another. The situation is another one when 
7i/2 < 7 1 r /2. Here, the eigenvalue trajectories avoid crossing in energy while the trajectories 
Ti and T 2 cross at certain values of a (Figs. [TJ e to h). 

The comparison of Fig. [TJ with complex u> = 0.05 (1 + i) and Fig. 2 in 9J with real 
u> = 0.05 shows the influence of lm(u). While the Tj(a) vary symmetrically around the 
critical value a cr when u is real, this is not the case when u is complex. In the last case, 
the difference Ti — T 2 blows up in the critical region. Furthermore, the Tj approach their 
asymptotic values in a relatively small parameter range of a when u is real in contrast to 
the case with complex u. These results will be discussed further in the following sections. 

In Fig. [2J we show the mixing coefficients \bij\ 2 of the eigenf unctions as a function of the 
parameter a corresponding to the eigenvalue figures shown in Fig. [TJ The mixing coefficient 
is complex, u = 0.05 (1 + 1). This figure can be compared with Fig. 5 in [9J calculated 
exactly with real 00 = 0.05. In both cases, the difference between the two curves \bij\ 2 is 
always 1 when 71/2 > 7^/2. This result corresponds to the fact that the two states will not 
be exchanged when 71/2 > 7^/2. The energy trajectories of the eigenstates cross freely and 
the widths of the two states are always different from one another, see Figs. [TJ a, b. 

At the exceptional point (where 71/2 = 7^/2), \by\ 2 — > 00 in both cases. When 71/2 < 
7f /2, the two curves |%| 2 coincide at a certain value of a. This happens for real (see Q]) 
as well as for complex (Fig. [2]) u, but is symmetrically only for real oj. At these points, the 
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FIG. 1: Energies Ei and widths Ti/2 (full lines) of N = 2 states beyond (a,b), at (c,d) and below (e to h) 
the exceptional point. The parameters of the subfigures are 71/2 = 1.2 (a,b); 7i/2 = 1.0 (c,d); 7i/2 = 0.9 
(e.f); 71/2 = 0.7 (g,h). Further parameters: e x = 1 - a/2; e 2 = a; 72/2 = I.I71/2; u> = (l+i) 0.05. The 
dashed lines show ei{a). 

two states are exchanged. 

In both cases 71/2 < 7f"/2 and 71/2 > 7f /2, the asymptotic values 1 and are reached 
in a smaller parameter range when u is real than in the case with complex u. Within this 
parameter range the two eigenfunctions are mixed. The range of mixing shrinks to one point 
when the two states really cross (at the exceptional point). It is especially large when two 
discrete states (with 7^ = 0) avoid crossing, as can be seen from Fig. 5 in [9] (where the 



11 




FIG. 2: The mixing coefficients \bij\ 2 of N = 2 states beyond (a), at (b) and below (c,d) the exceptional 
point. The parameters are the same as in Fig. [TJ 

calculations are performed exactly with real uj for all by). 

The parameter range in which two states are mixed due to the existence of an exceptional 
point in the neighborhood plays an important role in realistic systems. It will be discussed 
in the following sections with more than one avoided crossing for both real and complex u. 

IV. EIGENVALUES FOR N = 4 CROSSING LEVELS 

Let us first compare the calculations with N = 2 and N = 4 states and imaginary 
coupling u = 0.05 i between discrete states and environment of scattering states (Fig. 
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FIG. o: Energies Ei and widths r\/2 (full lines) of N = 2 states (a,b,g,h) and N = 4 states (c to f and i, j) coupled 
to K = 1 channel. The parameters of the subfigures are j t /2 = 0.5 (a to d); 7^/2 = 0.494, 0.498, 0.502, 0.506 (e,f); 
7l/2 = 0.45, 72/2 = 0.55 (g,h); 71/2 = 0.35; 72/2 = 0.45; 73/2 = 0.55; 74/2 = 0.65 (i,j). Further parameters: N = 2: ei = 
1 — a/2; e2 = a ; N = 4 : ei = 1 — a/2; e2 = 1.05 — a/2; es = 1.1 — a/2; e4 = a; ui = 0.05 i. The dashed lines show ei(a). 

In the calculations, the eigenvalue trajectories Ei(a) and Tj(a) are traced as a function of the 
parameter a which is varied in such a manner that one state crosses in energy one and three, 
respectively, other states. The results show two exceptional points in both cases appearing 
at the two different parameter values a™ and Figs. Oa, b for two states and Figs. 
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FIG. 4: The same as Fig. [3] but w = (1 + i) 0.05 (a,b,e,f); uj = 0.05 (c,d,g,h); N = 4; 7i = 0.5 (a to d); 
7,: = 0.35; 0.45; 0.55; 0.65 (e to h). 

EJc to f for four states. In between these two values, the widths bifurcate : the differences 
Ar/2 = \Ti/2 — r 2 /2| blow up. The two-level case is described analytically by equations 
( ITOl) to ( fl~6l) . The width bifurcation is completely reproduced in the numerical results. 

The parameter range \af — af\ is larger in the case with four states than in the other 
one due to the larger region in which avoided level crossings take place. Furthermore, the 
difference Ar/2 is larger in the case of 4 states than in the case of 2 states. In the case that 
all states have equal widths ji/2, the widths Tj/2 of only two states bifurcate also in the 
case with four states. The widths of the two other states remain unchanged (Figs. EJc, d). 
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FIG. 5: The same as Fig. EJbut N = 4 and K = 4; uj = (1 + i) 0.05 (a,b); w = 0.05 i (c,d); w = 0.05 
(e,f); 7 ,/2 = 0.5. 

This result holds also when the different r y%/2 differ slightly from one another (Figs. Ele, f). 
In this latter case, the two states avoid crossing at the two critical values af and af and 
the width bifurcation is almost the same as in the case with equal widths ji. 

The situation changes completely when the widths 7$ differ stronger from one another. 
When \ji/2 — ji±i/2\ > |e* — ej±i|, the states do not avoid crossing neither in the case with 
two states nor in the case with four states (Figs. Elg to j). The eigenvalue trajectories cross 
freely in energy and the widths do (almost) not bifurcate. The dynamics of the system is 
therefore completely different from that determined by the results shown in Figs. Oa to f. 

In Fig. UJ we show the eigenvalues Ei and Tj/2 for the case with four levels and one channel 
with complex and real coupling coefficients (u = 0.05 (1 + i) and oj = 0.05, respectively). 
Width bifurcation can be seen when uj is complex and the widths 7, are equal (Figs. HJa, 
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b), while the Tj = 7, are independent of a when u is real (Figs. Hlc, d). The results with 
different 7, differ from those obtained with imaginary coupling u = 0.05 i (Figs. |3jg to j) : 
the widths of two states cross at parameter values at which their energies avoid crossing 
(Figs. He to h). 

In Fig. [5J we show the results of calculations with equal number of states and channels, 
N = 4 and K = 4. Although usually K < N in realistic systems, the results allow us to 
receive a deeper understanding for the spectroscopic redistribution processes taking place in 
the critical region. The calculations are performed with the same parametric dependence of 
the energies as in the foregoing calculations with N = 4 and K — 1. Further assumptions: 
7i = 0.5 and co = 0.05 (1 + i), 0.05 i, and 0.05, respectively. 

The results show width bifurcation when u is complex or imaginary (Figs. [5Ja to d) 
which is however smaller than in the corresponding cases with K = 1 (Figs. HI a, b and Figs. 
EJc, d, respectively). Instead, the widths r, of the states are changed asymptotically due 
to the coupling to the different channels. When u is real (Figs. 0e,f), the widths Tj/2 are 
independent of a and are equal to ji/2 = 0.05 (i — 1 to 4). 



V. EIGENFUNCTIONS FOR N = 4 CROSSING LEVELS 

The coefficients bij defined in ( TTTj) determine the mixing of the eigenfunctions $j in 
relation to the basic wavef unctions <f)j. The mixing is nonvanishing in the neighborhood of 



the critical points of avoided level crossings as has been s 
for two crossing states and real coupling coefficients Uij 



rown in exact numerical calculations 
In the calculations of the present 



paper, the mixing of the wavefunctions is simulated by assuming a Gaussian distribution for 
the Uij, equation (}T8]) . It reproduces very well the exact results obtained in [9|. 

In the following, we consider the mixing coefficients for all wavefunctions $j the 
eigenvalue trajectories of which are studied in Figs. [3] to [5] as a function of the parameter a. 
In Fig. [61 the coefficients \bij\ 2 are shown for two and four states with imaginary coupling 
co = 0.05 % (see the corresponding eigenvalue trajectories in Fig. [3]). When the widths 7^ of 
all states are equal (Figs. [6ja, b), |&jj| 2 — > 00 at the two exceptional points af and af. In 
the parameter range af < a < af, the wavefunctions are mixed while ba — > 1 and b^ — > 
for i 7^ j beyond the two values af 2 - The states are completely (1:1) mixed in the range 
af < a < af in the two-level- case (Fig. Oa). The picture is more complicated in the four- 
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FIG. 6: The mixing coefficients | bi j | 2 of N = 2 and N = 4 states the eigenvalues of which are shown in Fig. [3] K = 1 
channel; £j = 0.05 i. 

level-case (Fig. EJb). Here, all states are involved in the redistribution taking place in the 
critical parameter range af < a < af. The exchange of the two states the widths of which 
remain unchanged according to Fig. |3jd can be seen from the energy eigenvalue trajectories, 
Fig. |3jc, as well from the mixing coefficients, Fig. |6jb. 
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FIG. 7: The mixing coefficients \bij\ 2 of N = 4 states the eigenvalues of which are shown in Fig. [4] K = 1 
channel; cj = (1 +i) 0.05 (a,c); w = 0.05 (b,d). 

The figures are similar when the widths ji of the states differ slightly from one another. 
Instead of exceptional points, there are avoided level crossings at af and af, see Fig. Elc 
for the four-level-case. The mixing coefficients show a dependence on a which is similar to 
that obtained for equal widths (Fig. Ob). 

When the widths 7« of the states i differ stronger from one another than in Fig. 0c, the 
eigenfunctions remain almost unmixed for all parameter values a (Figs. Eld, e). This result 
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FIG. 8: The mixing coefficients \bij\ 2 of N — 4 states the eigenvalues of which are shown in Fig. [5] K = 4 
channels; u = (l + i) 0.05 (a); cj = 0.05 i (b); u = 0.05 (c). 

corresponds to the almost constant width trajectories I\ as a function of a (Figs. |3jh, j). An 
exchange of states does not take place (Figs. Og, i). These results show the great influence 
of the exceptional points onto the dynamics of the system considered. 

The calculations in Figs. [3] and [6] are performed with imaginary coupling oj = 0.05 i. 
For comparison, we show in Fig. the mixing coefficients \hj\ 2 when u is complex and 
real, respectively, for the four-level-case (corresponding to the eigenvalue trajectories in Fig. 
HJ). When all 7^ are equal to one another and u is complex, the mixing coefficients \bij\ 2 
point to the participation of all states in the redistribution process in the critical parameter 
range (Fig. 0a). This corresponds fully to the corresponding eigenvalue trajectories that 
show level exchange as well as width bifurcation (Figs. HJa,b). When the coupling is real, 
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however, the wavefunctions of only two states are mixed while the wavefunctions of the other 
two states remain almost unchanged in the critical region (Fig. 0b). 

The mixing of the eigenf unctions is completely different from that discussed above when 
the widths 7$ differ more strongly from one another. A mixing of the wavefunctions appears 
only when the energy trajectories Ei(a) avoid crossing and the width trajectories Tj(a) cross, 
compare Figs. Hie to h with the corresponding Figs. [TJc, d. A width bifurcation does not 
take place. 

The mixing coefficients \bij\ 2 for the case with N = 4 states coupled to K = 4 channels 
are shown in Fig. |H] (parameters the same as in Fig. The differences between the 

cases with coupling to 4 channels to those with coupling to only 1 channel can be seen by 
comparing Fig. [SJa with Fig. 0a, Fig. IHJb with 0b and Fig. [HJc with Fig. [Tib. The mixing 
of the wavefunctions by coupling the system to four channels is more complicated than 
that by coupling to only one channel. This statement holds true also when oj is imaginary. 
Furthermore, an additional mixing is caused by the avoided level crossing at small a in the 
four-channel-case when u is complex. 

VI. EIGENVALUES FOR N = 10 CROSSING LEVELS 

We continue our studies with the matrix (JHJ) by choosing N = 10 states coupled to 1 
channel. The eigenvalue trajectories Ei(a) and Tj(a) are shown in Figs. [9] and QUI In 
these cases, one of the states crosses in energy the remaining nine states one by one. The 
nine states have the same parametric energy dependence of a and are shifted equidistantly 
relative to one another. 

In Fig. [9ja, b we show the results obtained with u = 0.07 i and 7$ = 0.5 for all states. 
The results are similar to those for the 4-level case (Figs. [3]c, d). Two exceptional points 
appear as well as width bifurcation : the width of one of the states is much larger while that 
of another one is much smaller than the widths of all the other states in the whole critical 
parameter range between the two exceptional points. 

In Figs. [9jc to f, the widths 7$ differ from one another. We compare the results obtained 
with an imaginary and with a real coupling constant, uj = 0.07 i and 0.07, respectively. 
When u is real (Figs. EJe, f), the energy trajectories Ei(a) are very regular: avoided level 
crossings appear one by one with all nine states. An exceptional point can be found only by 
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FIG. 9: Energies Ei and widths Ti/2 (full lines) of N = 10 states coupled to K — 1 channel. The parameters 
of the subfigures are jj2 = 0.5 (a,b); ji/2 = 0.50; 0.51; 0.52; 0.53; 0.54; 0.55; 0.56; 0.57; 0.58; 0.59 (c 
to f); u = 0.07 i (a to d); 0.07 (e,f). Further parameters: a = 1 - a/2; 1.1 - a/2; 1.2 - a/2; 1.3 - 
a/2; 0.9 -a/2; 0.8 -a/2; 0.7 -a/2; 0.6 - a/2; 0.5 - a/2; a. The dashed lines show e, (a). 

analytical continuation into the continuum (see The width trajectories I\(a) are less 
regular : they cross at several parameter values a. According to the results obtained for a 
smaller number of states in sections II III to [V} the widths do not bifurcate when u is real. 

When to is imaginary, the energy trajectories Ei(a) are also regular (Fig. 0c). Fig. [9ld 
shows nicely how widths continue to bifurcate such that the broad state remains broad and 
the narrow state remains narrow for all a in the critical region. Mostly, the widths of two 
states are different from one another when the energy trajectories cross. The states can, 
nevertheless, be exchanged as can be seen from Fig. [9jc. 

In Fig. [TUl we show results obtained with the imaginary coupling constant u = 0.05 i 
(being smaller than in Figs. [He, d) as well as those obtained with the larger values u = 
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FIG. 10: The same as Fig. Ubut lo = 0.05 i (a,b); uj = 0.08 i (c,d); u = 0.10 i (e,f); uj = 0.15 i (g,h). 

0.08 i, 0.10 z, 0.15 i (Figs. [10ja,b and Figs. [TOjc to h, respectively). We see width bifurcation 
similar as in Fig. [3jd for N = 4 states and in Fig. |9jb for 10 states with equal widths ji. In 
any case, width bifurcation appears in the whole critical parameter range between the two 
exceptional points. The width bifurcation occurring in the case of 10 levels is stronger than 
that with a smaller number of levels (compare section HV|) . The 9 states crossed by the state 
10 are exchanged in the critical parameter region as can be seen from the corresponding 
Figs. HOja, c, e, g. 

The results shown in Figs. [9] and [10] follow from calculations with a toy model being 
symmetrical in relation to the crossing states. In the critical parameter range, the states 
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FIG. 11: The same as Fig. [9] but e w = 0.15/(0.15 + a): uj = 0.07 i (a, b); w = (1 + i) 0.07 (c, d) and 
ui = 0.07 (e,f). 

are exchanged in all cases but 8 of them do not contribute actively to the width bifurcation 
when u is imaginary. This result follows from the high symmetry of the exceptional points in 
relation to the crossing states, i.e. from the linear dependence of the energy of the crossing 



state on the parameter a, eio(a) = a. The 8 states play the role of 'observers' simi 



arly as 



discussed in a toy model with 4 states (171] and in a realistic model with 3 states [22 ]. 

In order to exclude the high symmetry of the exceptional points in relation to nearby 
states we show in Fig. [11] the results of another version of the toy model. Instead of 
e io( a ) = a in Figs. 191 and [T0~| we use e w = 0.15/(0.15 + a) in Fig. [TTJ The width Tj of one 
of the states becomes large in the critical parameter range when u is imaginary or complex 
(Figs. [TTJb and d, respectively), while the widths of all the other states are much smaller. 
This result reflects the situation observed in many realistic cases (see e.g. the review js| and 



also [23]). A separation of a short-lived state from the other ones does not appear when u> 
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is real (Fig. Elf). 



VII. DISCUSSION OF THE RESULTS 

In sections II III to Wl\ we showed results obtained numerically for the eigenvalues Si = 
Ei — i/2 Ti and eigenfunctions $ ? by using the matrix (JSJ) (or (jSJ)) with 2, 4 and 10 levels, 
respectively. Only the energies e ? - of the states are varied as a function of a certain parameter 
a. The widths 7$ are assumed to be constant, the interaction between the levels z and 
j is either zero or simulated by the Gaussian distribution ( I18p where u is the same for all 
states (and also independent of a), and the self-energy terms are considered to be included 
into the diagonal matrix elements. Most interesting is the one-channel case for different N, 
e.g. the matrix ([6]) for N = 4. 

This toy model describes well the generic features of open quantum systems. The Wjj are 
complex. They stand for the interaction of the states via the environment (continuum of 
scattering wavefunctions) : Re(a>y) arises from the principal value integral (E]) and Im(o;jj) is 
the residuum (j3J). Both parts play an important role for the dynamics of the open quantum 
system at high level density where the resonance states overlap. The influence of Re(ujij) 
onto eigenvalues £3 and eigenfunctions $j is studied numerically exact for N = 2 levels in an 
earlier paper [9j . It causes the avoided crossing of the states which is well known for discrete 
as well as for narrow resonance states. 

Equations ( TT3|) to ( fl6|) show analytically the main difference between Re(cUij) and Im(wjj) 
in the 2-level case when the widths of the two states are equal, 71 = 72. According to ffT3]) . 
an exceptional point can be found only by analytical continuation into the continuum when 
Uij is real. Thus, the two states avoid crossing as it is very well known for discrete as 
well as for narrow resonance states [9(. In contrast, Im(u;y) causes two exceptional points 
according to equation (|14p . Most interesting is the width bifurcation arising in between the 
two exceptional points according to f fT6|) . 

The analytical results (fl3|) to ( fl6|) are well reproduced in our numerical calculations for 
the eigenvalues Si, Figs. |3ja, b. Also the eigenfunctions <3>j show the two exceptional points : 
\bij\ 2 — > 00 in approaching them, Fig. [6ja. In between the two exceptional points, the 
wavefunctions of the two states are strongly mixed. The mixing is 1:1 in the middle between 
the two singular points. Beyond the critical region between the two exceptional points, the 
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\bij\ 2 approach 1 when i = j and when i ^ j. Here, the two states may be exchanged, at 
most. The figures 3.a,b and 6. a represent the basic process of width bifurcation according 
to the analytical expressions (1151) and (1161) . 

In sections IIV[ IVl and I VI I of the present paper, numerical results for 4 and 10 states are 
shown under similar conditions, i.e. with equal or nearly equal 7, for all states : Figs. EJc 
to f and Figs. [6jb and c for the eigenvalues and eigenf unctions, respectively, of 4 states, 
and Figs. |9ja, b for the eigenvalues of 10 states. The development of width bifurcation as a 
function of increasing imaginary coupling vector u can be seen in Fig. [10] for iV = 10 states 
when the single widths 7^ are different from one another. The uniform width bifurcation in 
Fig. Ob appears at a smaller value of Im(u;) than in Fig. [TDJh due to the different values of 
the single 7^ in the last case. Generally, the width bifurcation is stronger when the number 
of states is larger. 

The influence of Re(uiij) onto the eigenvalue trajectories and eigenf unctions is also shown 
for iV = 4 and 10 (Figs. HJc, d for the eigenvalues of 4 states, Fig. Ob for the eigenf unctions 
of 4 states and Figs. EOe, f for the eigenvalues of 10 states). Under the influence of Re(cjjj), 
the states avoid crossing in energy and the widths do not bifurcate. When the calculation is 
performed with complex coupling vector (Re(uij 7^ 0) and lm(uij 7^ 0)), the widths bifurcate 
with some shift of the position of the maximum relative to that of the minimum (Figs. [Ha, 
b) and a larger parameter range of mixed wavefunctions (Fig. 0a). 

According to the results presented in sections [IV] to IVII the eigenvalue trajectories are 
strongly influenced by the values of the external (fixed) parameters when the resonances 
overlap and exceptional points determine the dynamics of the system. However, the eigen- 
value trajectories are almost independent of one another when the degree of resonance over- 
lapping is small. Examples are Figs. [3]g to j and [He to h for the eigenvalues of, respectively, 
two and four resonance states and the corresponding Figs. [6]d,e and [3c, d for the eigenfunc- 
tions of these states. These results show the strong influence of external parameters onto the 
dynamics of the system at high level density which is known from the study of realistic cases. 
An example is the enhanced transmission through microwave cavities when the formation of 
whispering gallery modes is supported by the manner the leads are attached to the cavity 
24]. 

Width bifurcation is directly related to the alignment of one of the N resonance states of 
the open quantum system to a decay channel (K = 1) with the consequence that it becomes 
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short-lived while other states become trapped (long-lived), i.e. more or less decoupled from 
the environment. Mathematically, the alignment of a resonance state is possible since the 
eigenf unctions of a non-Hermitian operator are biorthogonal and their phases (relative to 
those of the eigenfunctions of the other states) are not rigid in approaching an exceptional 
point, see equation (1121) . Width bifurcation appears therefore in our calculations when 
K = 1 and N = 4 or 10. This corresponds to realistic situations in which usually K < N . 
Width bifurcation does, however, not appear when N = 4 states are coupled to K — 4 
channels (Fig. [5]) although the wavefunctions are strongly mixed in the critical parameter 
region also in this case (Fig. [S]). 

An artifact of our model calculations is the assumption of equidistant energies of all 
but one state which are crossed by one state the energy of which depends linearly on the 
parameter a. Thus, the critical points are symmetrically in relation to the two neighboring 
states, and N — 2 states do not contribute actively to the width bifurcation taking place in 
the whole system. In Fig. [TTJ we show results obtained when the symmetry is somewhat 
disturbed (as it is usually the case in realistic systems). In such a case, the width of 
one state separates from those of all the other states. That means, all nine states remain 
trapped (decoupled) in the whole critical parameter range. This result corresponds to results 
obtained theoretically as well as experimentally in different realistic systems (see section 4 
in the review jsj). 

It should be added here, that asymptotically (beyond the critical parameter region) the 
states can be exchanged at most. This is the case, indeed, in all our calculations. All 
redistribution processes caused by the exceptional points take place only in the critical 
parameter region. This holds true for real as well as for imaginary coupling coefficients 
Uij and can be seen also in the mixing coefficients fry which approach 1 or when i = j 
and % 7^ j, respectively, beyond the critical parameter region. The two cases with real and 
imaginary u differ, however, by the length of the critical region. 

VIII. CONCLUSIONS 

The results of the present paper show that exceptional points cause width bifurcation 
because the coupling of the states via the environment (continuum of scattering wavefunc- 
tions) is complex. An example are the bound states in the continuum appearing at a finite 
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value of the coupling strength between system and environment when the interplay between 
internal and external interaction is taken into account, see section 4.4 of the review [8]. Of 
particular interest is the width bifurcation appearing at high level density. Here, many ex- 
ceptional points are near to one another. As shown in the present paper, width bifurcation 
causes, under this condition, a splitting of the system into two parts one of which exists in 
the short-time scale while the other one appears in the long-time scale. The two parts can 
not be observed together. The long-lived states occur as fluctuations in the short-time scale 
while the short-lived states perform a smooth background in the long-time scale. 

This phenomenon is known in literature from theoretical as well as experimental studies 
on different small quantum systems and systems equivalent to them (see examples in the 



review 



and in 



13, 
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26j). In many-body open quantum systems, it is called mostly 



dynamical phase transition. It is known also in optics where it is called superradiant phase 
transition according to Dicke 



mostly PT-symmetry breaking 



4j| . In PT-symmetric systems, the phase transition is called 



27]. 



Common to all these studies is that the dynamical phase transition is very robust when 
the necessary conditions are fulfilled, i.e. when the level density is high and the dynamics of 
the system is determined by exceptional points. Although the existence of exceptional points 
is decisive (as shown in the present paper), the dynamical phase transition does not appear 
at the parameter value which gives the position of the exceptional point itself. Instead, the 
phase transition occurs in the neighborhood of one (or several) exceptional points. Such 
a result is known from different experimental studies on concrete realistic systems. It is 
discussed in e.g. 28j. The results of the present paper show that the critical parameter 
range is determined by the distance between (at least) two exceptional points. 

The dynamical phase transition appearing in small quantum systems at high level density, 
causes a new understanding of time as discussed in 11] . Time which is characteristic of the 
system, is inverse proportional to the widths Tj and related therefore directly to the non- 
Hermitian part of the Hamilton operator. The widths do not increase limitless as shown in 
the present paper. Instead, the system is dynamically stabilized : the lifetimes of the states 
of the system are increased and the system as a whole is stabilized by ejecting the short-lived 
state from the system at the dynamical phase transition. Thus, the time characteristic of 
the system is bounded from below in a similar manner as the energy. Beyond the dynamical 
phase transition, all states have lost their original spectroscopic features and the number of 
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states is reduced. 

This fact is very well known in nuclear physics (although very seldom interpreted in this 
manner). The properties of compound nuclei are described well by the Unified theory of 
nuclear reactions developed by Feshbach which contains both the short-lived direct reac- 
tion part and the long-lived compound nuclear reaction part. However, Feshbach introduces 
statistical assumptions in order to describe the long-lived compound nuclear states. This is 
in contrast to the calculations in the present paper where the long-lived states are described 
without any statistical assumptions. Instead the long-lived states are shown by us to arise 
from width bifurcation causing a dynamical stabilization of the system at high level density, 
a phenomenon that is caused by exceptional points. 

Many of these results seem to be counterintuitive. They allow us, however, to explain 
some unexpected experimental results. For example, the crossover from the mesoscopic 
to a universal phase for transmission in quantum dots observed experimentally [3] can be 
explained qualitatively by the formation of the short-lived resonance state at the dynamical 
phase transition. The results of the present paper support this interpretation [12j of the 
experimental results. Another example is the experimental observation j^] of non-statistical 
effects in nuclear reactions on middle-heavy nuclei. The data show directly the formation 
of different time scales in the system at high level density. Other examples are discussed in 

The long-lived states beyond the dynamical phase transition are strongly mixed and show 
chaotic features. The random matrix theory is surely applicable only after averaging over 
different decay channels (as usually done), i.e. it is not applicable to the description of small 
systems coupled to one channel (and, respectively, to two channels in the case of trans- 
mission through the system). A theoretical analysis of the spectra (without any statistical 
assumptions) by restricting to the results obtained from one decay channel according to the 
recent experimental data on compound nuclei [5| is not performed up to now. On the basis 
of the results of the present paper, it can be stated today only the following. The mix- 
ing of the long-lived states (beyond the dynamical phase transition) is caused by complex 
many-body forces via the continuum of scattering wavefunctions (simulated by the Uij in 
our calculations) and not by two-body forces. This may provide an explanation of the fact 
why compound nucleus spectra (after averaging over different channels and beyond decay 
thresholds) may be described by a Gaussian orthogonal ensemble (containing many-body 
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forces), but not by a two-body random ensemble. This point should be studied in future in 
more detail. 

In order to prove the role of exceptional points in an open many-level quantum system it is 
highly interesting to study experimentally time symmetry breaking caused by the influence of 
a nearby state onto an exceptional point. Theoretical studies with a symmetrical influence 
(as in most calculations of the present paper) can surely not describe the properties of 
realistic systems with broken time symmetry. Symmetry breaking influences not only width 
bifurcation (as shown in Fig. [TTj) but will prove, above all, the irreversible processes caused 
by exceptional points in open quantum systems. These processes are decisive, among others, 
for the dynamical stabilization of quantum systems and the formation of quantum chaos. 

Finally we remark that the sensitive dependence of the dynamics of an open quantum 
system on the value of external parameters (as shown in the present paper) can be used 
in order to construct small systems with desired properties. This point is important for 
applications. 

* email: eleuchh@iro.umontreal.ca 
** email: rotter@pks.mpg.de 
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